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Globular clusters are highly efficient radio pulsar factories. These pulsars can be used as precision 
probes of the clusters’ structure, gas content, magnetic field, and formation history; some of them 
are also highly interesting in their own right because they probe exotic stellar evolution scenarios 
as well as the physics of dense matter, accretion, and gravity. Deep searches with SKA 1-MID and 
SKA 1-LOW will plausibly double to triple the known population. Such searches will only require 
one to a few tied-array beams, and can be done during early commissioning of the telescope — 
before an all-sky pulsar survey using hundreds to thousands of tied-array beams is feasible. With 
SKA2 it will be possible to observe most of the active radio pulsars within a large fraction of 
the Galactic globular clusters, an estimated population of 600 — 3700 observable pulsars (those 
beamed towards us). This rivals the total population of millisecond pulsars that can be found in 
the Galactic field; fully characterizing it will provide the best-possible physical laboratories as 
well as a rich dynamical history of the Galactic globular cluster system. 


Advancing Astrophysics with the Square Kilometre Array 
June 8-13, 2014 
Giardini Naxos, Italy 


* Speaker. 


© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 








Pulsars in GCs 


J. W. T. Hessels 


1. Introduction: The Science of Globular Cluster Pulsars 

Compared with the Galactic held, globular clusters (GCs) contain large numbers of millisec¬ 
ond pulsars (MSPs) per unit stellar mass — for recent reviews see: Camilo & Rasio (2005); Ran¬ 
som (2008); Freire (2013). The stellar densities in the cores of GCs are orders-of-magnitude larger 
than in the Galactic plane; these conditions create a unique environment for stellar collisions and 
interactions (e.g. Sigurdsson & Phinney 1993, 1995). As a result, neutron stars can swap in new 
companion stars and form exotic binary (and trinary) systems (e.g. Thorsett et al. 1999; Sigurdsson 
et al. 2003), in which the neutron star is often spun-up to become an MSP through a period of 
mass and angular momentum transfer (Alpar et al. 1982; Radhakrishnan & Srinivasan 1982). Fur¬ 
thermore, GCs have half-light radii of only a few arc-minutes or less, making them easy targets for 
directed searches. Deep searches of the Galactic GC system have thus far discovered 144 pulsars in 
28 clusters^ This includes many remarkable systems, e.g.: the fastest-spinning neutron star known 
(Hessels et al. 2006), potentially the most massive neutron star known (Freire et al. 2008a,b), a 
double neutron star system (Prince et al. 1991), a triple system with a white dwarf and Jupiter-mass 
companion (Thorsett et al. 1999; Sigurdsson et al. 2003), as well as many highly eccentric binaries 
(Ransom et al. 2005; Freire et al. 2004, 2008a; Lynch et al. 2012) and systems with ‘odd’ orbits or 
stellar companions (Lyne et al. 1993; D’Amico et al. 2001). 

Population simulations, based on the results of deep pulsar searches in the past decade, show 
that the Galactic GC system harbors many more pulsars still to be discovered (Bagchi et al. 2011; 
Chennamangalam et al. 2013; Turk & Lorimer 2013). In particular, once extrapolated to the total 
sample of 157 known GCs in the Milky Way, the most recent study (Turk & Lorimer 2013) pre¬ 
dicts a population range of potentially observable pulsars (i.e. those beamed towards us) between 
600 — 3700 (95% confidence level). The distribution of this population among the various known 
sub-classes of pulsars (isolated, binaries, trinaries, eclipsing, etc.) has not been carefully explored 
(though see Verbunt & Freire 2014). Nonetheless, past experience indicates that any significant 
growth in the known population inevitably leads to the discovery of some exotic and rare new 
kinds of pulsar systems (e.g. the pulsar triple system PSR J0337-I-1715; Ransom et al. 2014). In 
this respect, pulsars in GCs are of special interest because the peculiar environment of the core 
of a GC provides the conditions for producing systems that are unlike anything that can be easily 
formed in the Galactic field (e.g. Freire et al. 2004), thus opening the possibility of detecting ‘Holy 
Grail’ systems, like the first-ever MSP -i- black hole binary or even an MSP + MSP binary. The 
known GC pulsars, and those that the SKA will discover, can be studied very efficiently, because 
they are all located close to each other on the sky. For example, a single observation using no more 
than a half dozen tied-array beams could time all 34 known pulsars in Terzan 5. 

The science of GC pulsars falls into two main themes: 

i) Individual systems that can be used for investigations in fundamental physics, like the 
equation-of-state of neutron-rich matter at super-nuclear densities and/or tests of gravita¬ 
tional theories (Ransom et al. 2005; Freire et al. 2008a,b; Demorest et al. 2010; Kramer et al. 
2006) — see also other relevant chapters in these proceedings (Keane et al. 2014; Shao et al. 2014; 

^See http://www.naic.edu/~pfreire/GCpsr.html for an up to date catalog. 
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Watts et al. 2014). Such systems can also provide novel insights for understanding stellar evolution 
and accretion physics (Archibald et al. 2009; Papitto et al. 2013; Bassa et al. 2014; Stappers et al. 
2014). Although these science cases can also be addressed using pulsars in the Galactic held, the 
radically different environment of a GC favors the formation of exotic objects, the peculiar param¬ 
eters of which may be even better suited for certain applications or constraining physical models. 

Ultra-Fast Rotators: GCs harbor some of the fastest-spinning neutron stars known, including 
the 716-Hz record holder in Terzan 5 (Hessels et al. 2006). Though such sources are rare (Hessels 
et al. 2007), doubling or tripling the known population provides great prospects for pushing towards 
even faster rotation rates — perhaps even a sw/^-millisecond pulsar, which could strongly constrain 
the maximum neutron star radius and hence the equation-of-state (Hessels et al. 2006; Lattimer & 
Prakash 2001). Unlike in the Galactic held, an MSP in a GC can in principle experience multiple 
episodes of recycling (due to dynamical encounters) and hence could be more effectively spun up 
to its limiting rotational period (as well as growing in mass). This remains to be shown directly, 
however. 

Neutron Star Masses: Many (at least 23) highly eccentric {e > 0.01) pulsar binary systems exist 
in GCs. For these, the neutron star mass can be constrained with the help of periastron precession 
(e.g. Freire et al. 2008a). In the Galactic held, MSP orbits are almost exclusively extremely circular 
{e < 10“"^), and prohibit such measurements in most cases. Besides the obvious implications of the 
highest-mass sources for constraining the neutron star equation-of-state (Antoniadis et al. 2013; 
Demorest et al. 2010), mapping the full MSP mass distribution is an important probe of their 
formation in supernovae, and their later ‘recycling’ through accretion. Once enough NS masses 
have been measured, the maximum NS mass should become apparent. 

Evolution and Accretion: Some individual systems also pose interesting evolutionary puzzles, 
which touch on accretion physics, accretion-induced collapse, and stellar exchanges. For example, 
a handful of slow, apparently young pulsars are also known in GCs (Lyne et al. 1996; Boyles et al. 
2011), providing interesting cases for exploring alternatives to the typical core-collapse supernova 
channel for forming neutron stars. The hrst MSP with a bloated main-sequence companion was 
found in a GC (D’Amico et al. 2001). This is now the prototype for a known population of a dozen 
such sources, some of which are seen to switch between radio MSP and low-mass X-ray binaries in 
quiescence and/or outburst (Archibald et al. 2009; Papitto et al. 2013; Bassa et al. 2014; Stappers 
et al. 2014). In fact, there are now 21 eclipsing pulsars (so-called ‘black widows’ and ‘redbacks’) 
in GCs, meaning that there are excellent prospects for studying more such transitions in the future. 

Multi-wavelength Studies: Deep multi-wavelength observations of GCs in X-rays (primarily 
using Chandra-, see e.g. Grindlay et al. 2001; Pooley et al. 2002; Grindlay et al. 2002; Heinke et al. 
2003a,b, 2005, 2006; Bogdanov et al. 2006; Eisner et al. 2008; Bogdanov et al. 2011), y-rays (pri¬ 
marily using Fermi-, see e.g. Freire et al. 2011), and optical (primarily using HST, see e.g. Edmonds 
et al. 2001; Bassa et al. 2003, 2004; Pallanca et al. 2010, 2013, 2014) provide complementary in¬ 
formation on either the pulsar’s magneto spheric emission, intra-binary emission (perhaps from a 
shock), or the companion itself. The next generation of extremely large telescopes (EETs) may pro¬ 
vide large numbers of radial velocities for constraining the mass ratios of the stellar components, 
as in Cocozza et al. (2006). Conversely, identifying high-energy GC sources as radio pulsars helps 
better understand the zoo of objects that can be created (e.g. cataclysmic variables, low-mass X-ray 
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binaries, etc.) and hence the stellar evolution history of the cluster itself (Heggie & Giersz 2008). 
For GC pulsars we have an independent measure of the distance (and often the reddening), which 
can allow stronger constraints on several pulsar+companion quantities measured in optical. X-rays 
or 7 -rays. The Fermi satellite is revealing 7 -ray emission from several GCs (e.g. Abdo et al. 2010). 
Comparing the total 7 -ray emission with the number and properties of the known radio MSPs in 
the largest possible sample of GCs will be a new tool for investigating the as yet poorly understood 
high-energy emission mechanisms of the MSPs (e.g. Harding et al. 2005), as well as constraining 
the MSP radio beaming factor (i.e. the fraction of the sky swept by the MSP radio beam). Fermi is 
not only detecting several GCs, but has started detecting individual pulsars, like PSR B1820—30A 
(Freire et al. 2011) and PSR B1821—24A (Johnson et al. 2013). 

ii) Using the cumulative population of pulsars to probe the structure, proper motion, dynami¬ 
cal status, magnetic field, and gas content of the cluster itself (Phinney 1993; Meylan & Heggie 
1997; Pooley et al. 2003). Pulsar-determined quantities can constrain the still largely obscure GC 
evolution, as well as its relation with the history of the Galaxy. For example, the pulsar-abundant 
cluster Terzan 5 contains at least two stellar populations of differing ages, and is possibly the pris¬ 
tine remnant of a building block of the Galactic bulge (Ferraro et al. 2009). The same advances 
in computing that will allow the creation of the SKA will also allow advances in numerical sim¬ 
ulations of GCs. Simulations of the stellar and dynamical evolution of stars in a GC will predict 
the observable population of neutron stars in each cluster, and this can be directly tested against 
observations. Neutron stars, as some of the heaviest objects in GCs, play a vital role by acting as 
an energy reservoir to counter the gravitational collapse through the formation and disruption of 
binaries. The SKA will provide a census by looking for those neutron stars seen as pulsars, and 
their spin and binary properties will give important clues to the dynamical state of the cluster (e.g. 
Verbunt & Freire 2014). 

Cluster Potential: Acceleration in the cluster potential affects the pulsars’ spin derivatives 
and binary orbital period derivatives, which can be used to probe the cluster potential in a very 
direct and unique way, possibly revealing the presence of an intermediate mass black hole at the 
core. Strong constraints can also be set on the inner mass-to-light ratio by finding a very centrally 
located pulsar (D’Amico et al. 2002). Thus far, there is no certified intermediate-mass black hole 
(IMBH) known to exist, anywhere in the Universe. Finding such an object would lead to a wealth 
of studies to understand its formation and characteristics, e.g. testing if the velocity dispersion of 
nearby stars in the GC is related to the black hole mass, as is seen in galaxies over a few orders- 
of-magnitude in mass. The discovery of an MSP orbiting an IMBH would give a unique chance to 
directly measure the black hole spin (Liu et al. 2014, see also Liu 2012, PhD, U. of Manchester). 

Proper Motion: The ensemble of pulsar proper motions, measured through pulsar timing, can 
be used to determine the cluster’s proper motion and hence to infer its orbit in the Galactic gravi¬ 
tational potential (to some extent optical studies already measure GC proper motions, although not 
very precisely; though this will change with GAIA). While radio can provide accurate proper mo¬ 
tion, optical instruments (especially 30-meter-class telescopes) can provide accurate radial motion. 
The combination produces a 3D velocity vector for the GC. When this is done for a sizable num¬ 
ber of GCs, it will be possible to constrain the Galaxy’s gravitational potential well. For nearby 
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clusters, there is also the chance to detect the peculiar motions of the MSPs within each cluster. 
If they are in binary systems, then also the optical radial motion and hence the 3D motion of the 
population of binary MSPs in a GC can be determined. Again this will be a handle for constraining 
the GC potential well, especially in the core, where the relatively massive MSPs (compared to the 
mean stellar mass) typically reside. 

Intra-cluster Medium: Differential dispersion and rotation measures (DMs and RMs) between 
the detected pulsars can be used as highly sensitive, unique probes of the intra-cluster ionized 
medium and magnetic held (Freire et al. 2001b). This can teach us about the stellar winds of old 
stars releasing plasma in the cluster and/or the interactions between this plasma and the winds of 
MSPs and hot stars. Furthermore, the existence of this plasma, together with sensitive X-ray limits 
on accretion luminosity from a central source, can be used to place an upper limit on the mass of a 
possible central black hole. As to the putative magnetic helds in GCs, these will be an additional 
ingredient in attacking the difficult problem of justifying the occurrence of large-scale magnetic 
helds in the Universe. 

Interaction History: The types of pulsars that are found — e.g. their spin-period distribu¬ 
tion, locations in the cluster, and the fraction of binary versus isolated pulsars — can vary quite 
drastically from cluster to cluster, and encodes information about the dynamical history of the GC 
(Verbunt & Freire 2014). Here we will learn about the still largely unknown stages of a GC’s 
evolution, with particular emphasis on the process leading to core-collapse. There are indications 
that the situation is quite different with respect to what was believed until recently; for the relevant 
theory and observations see Fregeau (2008) and Pooley (2010), respectively. 

There are 157 known Galactic GCs (Harris 1996)^, most of which have now been searched for 
pulsars. Since the earliest handful of GC pulsar discoveries (Lyne et al. 1987, 1988; Manchester 
et al. 1990; Anderson et al. 1990; Manchester et al. 1991; Kulkarni et al. 1991), the known popula¬ 
tion has exploded in the last decade. This has been thanks to a number of large surveys with Parkes 
(e.g. Camilo et al. 2000; Ransom et al. 2001; Possenti et al. 2003), Arecibo (e.g. Hessels et al. 
2007), and the Green Bank Telescope (GBT; e.g. Ransom et al. 2004, 2005). All of these surveys 
greatly benehtted from new, wide-band pulsar data recorders that provided big leaps in time and 
frequency resolution, critical for identifying the fastest rotating pulsars. The hrst major crops of 
GC pulsars were in M15 (Anderson 1993) and 47 Tuc (Manchester et al. 1990, 1991; Camilo et al. 
2000; Freire et al. 2003, 2001a). Deep synthesis imaging underlined the potential for discovery 
(Fruchter & Goss 2000). For example, deep imaging of Terzan 5 (Fruchter & Goss 2000) spawned 
great interest in that cluster, which later resulted in a record 34 known MSPs being found (Ransom 
et al. 2005; Hessels et al. 2006). The discovery rate has recently fallen off, as we have reached the 
limits of what the current generation of large single dishes can do. As we now discuss, SKAl can 
provide the next boom in GC pulsar discoveries, and this can be achieved in the very early days of 
the telescope, even with a subset of the total collecting area. 

2. Observing Pulsars in Globular Clusters with SKAl 

The search sample: SKAl will provide unprecedented sensitivity for targeted GC searches 
^Last revision from December 2010 available at http://physwww.mcmaster.ca/~harris/mwgc.dat. 
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in the southern skies, which is fortunately also where the majority of the most massive, densest, 
and hence most pulsar-rich GCs can be found, i.e. in the Galactic bulge. Compared with the 
greater observing and processing challenge of an all-sky search, which will require the formation 
of hundreds to thousands of tied-array beams, deep targeted searches of these limited helds-of- 
views will provide early commissioning science by quickly hnding dozens, and eventually a few 
hundred pulsars in these clusters (once sufficient collecting area is in place, even if only 50% of 
the total, and at least a basic tied-array mode is available; we also note that even an incoherent sum 
of > 50 dishes would already provide higher sensitivity for those clusters that are currently only 
visible to Parkes). Follow-up timing observations will also be facilitated by the small required field- 
of-view. We know that we are currently only sampling the tip of the pulsar luminosity distribution 
in these predominantly distant (> 4kpc) stellar systems. In fact, recent investigations (Bagchi et al. 
2011; Chennamangalam et al. 2013) conhrmed that the luminosities (dehned as Ly = Sd^, where 
S is the observed flux density at central frequency v in MHz, and d is the GC distance) of the 
pulsars observed in a GC can be reproduced as the bright tail of either a log-normal distribution 
— with parameters compatible with the luminosity functions (LFs) inferred for the pulsars in the 
Galactic held — or a power-law distribution with index ^ — 1 (in agreement with earlier results by 
McConnell et al. 2004; Hessels et al. 2007), with the former functional form providing a slightly 
better match to the available data. For both the assumed LFs, in a large range of not too weak 
luminosities (typically above 0.5 mJy kpc^), the flux density distribution follows dlogN/dlogS^ 
0.5 — 1 for any given cluster. That implies that a large increase in sensitivity will automatically 
bring many new sources within the reach of detection and will also provide higher-precision studies 
of the (mostly very faint, 5i4oo ^ 20 /iJy) sources that are currently known. This will also enable, 
e.g., more mass measurements and equation-of-state constraints using already known sources. As 
we outline below, SKAl can plausibly increase the known population of 144 pulsars by a factor 
of two to three and, most importantly, excellent laboratories for studying dense matter physics and 
gravity will be discovered. In some cases, it may be possible to detect all the active radio pulsars 
in a given cluster, providing a unique view of the star formation history and interactions over the 
cluster’s lifetime. In most cases, however, the SKA2 will be needed to reach the full observable 
sample. 

Based on population synthesis simulations (Turk & Lorimer 2013), a conservative total of 
600 — 3700 detectable (i.e. beamed towards us) pulsars in the Galactic GC system can be inferred. 
This rivals the total expected population of MSPs in the Galactic field, estimated at a few tens of 
thousands of sources (Faucher-Giguere & Kaspi 2006; Lorimer 2008). Importantly, while Galactic 
field pulsars are spread across the full 41,000 sq. degrees of the celestial sphere, the total area 
required to search all known Galactic GCs combined constitutes only ~ 1 sq. degree! Of the 157 
known Galactic GCs, there are 154/152 GCs visible to SKAl-MID/SKA 1-LOW for at least 2hrs, 
assuming elevation limits of 15/30 deg, respectively. For 8-hr integrations, there are still 141/113 
clusters available for SKAl-MID/SKA 1-LOW, respectively. Of these, 28 have at least one known 
pulsar, and hence the DM to the cluster is also known (a fact that makes searching significantly 
easier). 

MSPs are in general weak radio sources (phase-averaged flux densities of 5*1400 ^ 1 mJy), and 
GCs are typically at > 4kpc. Raw sensitivity and long dwell time is thus crucial for finding the 
weakest cluster pulsars, though long dwell time also adds complications, as we discuss below. The 


6 



Pulsars in GCs 


J. W. T. Hessels 


weakest known GC pulsars have *Si 4 oo ^10 /iJy and were discovered in multi-hour integrations 
using the GBT and Arecibo with several hundred MHz of bandwidth at 1.4 or 2 GHz. This sets the 
bar that SKAl must surpass. 

Searching for binaries: Greater instantaneous sensitivity is not only important for hnding 
weaker sources, but also because more than two-thirds of the known sources are in a binary system. 
Orbital motion smears the pulsar signal over multiple Fourier bins, and can make such sources 
undetectable unless this is corrected for (Ransom 2001). For this reason, it is likely that the intrinsic 
fraction of binary MSPs in clusters is larger than the observed fraction. Typically, binary pulsars are 
discovered using an ‘acceleration search’, which approximates the Doppler shift of the signal as a 
constant drift in the frequency domain (Ransom 2001). This approach is relatively computationally 
efficient; full searches of all Keplerian orbital parameters (even in the case of a circular orbit, where 
only three orbital parameters are needed) are currently not tractable for general use. Acceleration 
searches are only valid when the integration time of the full observation, Fobs, is less than 10% of the 
orbital period, Porb- In other words, one cannot necessarily gain sensitivity to binary pulsar systems 
simply by integrating for a longer time (though some analytical techniques do exist; see Ransom 
2001). The shortest known orbital period of a GC pulsar is Porb = 1.5 hr. Such a system can thus be 
found in a linear acceleration search of a 10-min data set. In the following, we therefore use 10 min 
as a fiducial integration time for searches for binary pulsars in compact orbits, while we use 2 hr as 
an appropriate integration time for longer-period binaries or isolated pulsars. The assumption of a 
constant acceleration is useful even for eccentric systems, although searches over a ‘jerk’ parameter 
may be helpful (Bagchi et al. 2013). With the constant-acceleration assumption, highly eccentric 
systems are in fact more easily detectable at most orbital phases than low-eccentricity systems^. 

Propagation effects: As in all pulsar surveys, interstellar propagation effects can strongly 
limit detectability, especially at the shortest spin periods, where residual smearing in time due to 
uncorrected dispersive or scattering delay can broaden the pulse in time to the point where it is 
no longer detectable. Many of the densest, most massive GCs known are located in the Galactic 
bulge. The line-of-sight DM is often large (> 200pccm“^), as is the expected scattering. To 
mitigate these effects, one can observe at higher radio frequency. Taking the competing effect 
of the typically steep spectra of pulsars {S oc where 1 < a < 3; Maron et al. 2000; Bates 
et al. 2013) into account, it turns out that the 1.4 — 2.0 GHz band is well suited for searching GCs 
with DMs greater than ~ 100pccm“^. Roughly 100 of the SKA-visible clusters have expected 
DM> 100pccm“^ (according to the NE2001 model of Cordes & Lazio 2002), and are excellent 
targets for SKAl-MID, which provides the maximum instantaneous sensitivity in the 1.4 — 2.0 GHz 
band. For the remaining ~ 60 SKA-visible clusters with lower expected DMs, SKAl-LOW (and/or 
SKAl-MID at 800 MHz in some cases) presents an exciting opportunity, as we outline below. 

Array configuration: GC pulsar searches require the SKA’s tied-array observing mode be¬ 
cause the data must be recorded with at least ~ 50 /is time resolution. Maximum possible sensi¬ 
tivity can be achieved if the data are coherently dedispersed at the average cluster DM (in the case 
of GCs with known pulsars) or at least at a few trial DMs for clusters with no known sources. In 
practice, this will likely require the VLBI-like ability to record raw voltages from the tied-array. 
An important consideration is the sky area that must be covered for such searches. The known 

^See Madsen 2013, MSc, UBC; https://circle.ubc.ca/handle/2429/44897 
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Galactic GCs typically have core radii of a few arc-seconds to one or two arc-minutes and half- 
light radii of ~ 0.5 — 3 arc-minutes (Harris 1996). GC pulsars congregate to a large extent towards 
the cluster core, and almost all of the known GC pulsars are within one arc-minute of the cluster’s 
optical center-of-light. The 400-m (radius) core of SKA 1-MID will provide a tied-array beam with 
FWHM of 0.9 arc-minute at 2 GHz, meaning that only a single beam will be needed to detect and 
discover the majority of the GC pulsars using this setup. In order to catch the minority of pulsars 
further from the core, which are still very interesting for understanding the dynamical history of 
the cluster, a more than sufficient coverage can be achieved with hexagonally packed grids of 7 
or 13 beams. For the 600-m (radius) core of SKA 1-LOW, a single 3 arc-minute tied-array beam 
at 350 MHz will be sufficient. Thus, searching for GC pulsars is an excellent early science and 
commissioning project, which will naturally lead the way to an all-sky pulsar survey requiring 
hundreds to thousands of tied-array beams. More ambitious searches using a larger fraction of the 
SKA 1-MID/LOW antennas will be more challenging. For example, to double the instantaneous 
sensitivity provided by SKA 1-MID requires going to a radius of 4000 m from the array center, and 
thus 100 times more beams to cover the same area of sky. This approach is likely more interesting 
for follow-up than for initial searches. For example, once new pulsars are found and sufficiently 
well localized, the further timing observations that are needed to extract the science (e.g. precision 
astrometry, orbital elements, proper motion, post-Keplerian effects, etc.) can be done by using one 
ndiUosN, full-array tied-array beam on each source. 

Achievable sensitivity: The 400-m (radius) core of SKAl-MID, operating from 1.4 — 2.0 GHz 
can achieve an rms noise of 4.0/1.2 /iJy for lO-min/2-hr integrations. Using a digital beam-former 
and back-end similar to those sketched in the SKA Baseline Design (i.e. capable of keeping the 
smearing effects of interstellar dispersion smaller than the adopted sampling time of 50 /is for all 
DMs < 1000pccm“^), that rms noise translates into a limiting sensitivity (at S/N = 8) of ^ 12/4 
/iJy for a recycled pulsar spinning at a ^ 1 ms period and having an intrinsic pulse duty cycle 
of 10%. Those sensitivity limits halve for the GCs where 8hr-long tracks with SKAl-MID are 
possible (this fortunately includes some of the most promising targets)"^. SKAl-MID can reach a 
factor of ~ 3 deeper than the state-of-the-art searches performed using the GBT and the GUPPI 
back-end (DuPlain et al. 2008). Only 32 clusters within declination, 5, of —2deg< 5 < +38 deg 
are visible to Arecibo. For most of the 30 clusters south of 5 = —40 deg the best limits have 
been reached by Parkes; SKAl-MID can go ~ 11 times deeper in these cases. The 600-m (radius) 
core of SKAl-LOW, operating from 250 — 450 MHz, can achieve ^ 100/30 /iJy sensitivity for 10- 
min/2-hr integrations (depending on the line-of-sight; for clusters in the Galactic bulge, Tsky will be 
significant). For comparison, the all-sky GBNCC survey at 350 MHz achieves ~ 1 mJy sensitivity 
(Stovall et al. 2014). 

Available discovery space: Figure 1 illustrates the discovery space that is open to SKAl-MID 
for these kind of pulsar searches (a similar analysis can be made for SKAl-LOW). We define the 
quantity F if to represent the growth of the probed area of the pulsar luminosity function (LF) of 
a given GC. This assumes a survey performed with SKAl-MID at 1700MHz central frequency, 

is important to note that longer integrations do not necessarily improve sensitivity to binary pulsars because only 
a coherent search of the possible orbital parameters can fully recover the Doppler-smeared signal. Nonetheless, by the 
time that the SKA is operational, it may be possible to use more advanced binary search techniques than those typically 
used at present. 
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SKAl-Mid - Factor of improvement r,^p in the search in GCs 



Figure 1: A histogram of the distribution of Flf, which represents the growth of the probed area of the 
pulsar luminosity function, for all GCs visible to the SKA. Here we assume use of SKAl-MID in 10-min, 
2-hr, and 8-hr observations, compared to state-of-the-art observations with Arecibo, GET, and Parkes. See 
text for further details. 


with 600 MHz bandwidth and other survey parameters like those in the SKA Baseline Design. In 
particular, the histogram of Fig. 1 is obtained by comparing the sensitivity of this SKA 1-MID survey 
with that of the best possible GC searches already carried out or still ongoing. For each cluster, 
we assume that the best-possible search can be conducted by either Arecibo, Parkes, or GET, 
depending on the declination of the cluster (Arecibo is the best telescope for all clusters within 
—2deg< d < +38 deg, Parkes is the best telescope for all clusters 5 < —40 deg, and the GET is 
the best available for 5 > —40 deg and outside the Arecibo declination range). For the purpose 
of comparison, all flux density limits have been scaled to 1400 MHz using a pulsar spectral index 
a = —1.7. We also assume, for simplicity, that the back-ends used in these various surveys are 
equally as good at minimizing the effects of dispersion. Finally, also the effects of scattering in the 
ISM have been assumed to be the same at the central frequencies of SKA 1-MID and at the reference 
frequencies of the other telescopes 1400 MHz for Parkes and Arecibo, and ~ 2000 MHz for 
GET). Using these assumptions, and an effective pulsar duty cycle of 25% (which is a typical for 
millisecond pulsars), we calculated survey flux density limits. The next step in producing Fig. 1 was 
to feed the calculated sensitivity limits to a log-normal pulsar LF with mean (in units of mJy kpc^ 
expressed in a logarithmic base-10 scale) of —1.1 and dispersion 0.9, which is known to reproduce 
the observed data (Bagchi et al. 2011; Chennamangalam et al. 2013). 

We considered the distributions of F^f for 3 hypothetical surveys looking at each GC for 10 
min, 2 hr, and 8 hr (color-coded in red, blue and green), respectively. It is evident that, for most 
GCs, SKA 1-MID provides up to a factor 5 improvement in probing the LF. The average Tip over 


9 
































Pulsars in GCs 


J. W. T. Hessels 


the entire population is between 3 and 4^. The higher average Y^p for shorter integration times 
is due to the shape of the log-normal LF, which is steeper at higher luminosity. In other words, a 
given increment in sensitivity produces higher values of F pp when only the brightest part of the LF 
is sampled, which typically occurs for short pointings. This effect also largely explains the second 
peak at Tpp ^ 10, which results from the ^ 10 very distant GCs that are visible from Parkes. The 
increment in sensitivity of SKA 1-MID with respect to Parkes is large and, moreover, if a pulsar 
should be discovered in those GCs by a shallow survey, its luminosity will have to be very high 
and thus Tpp is correspondingly very high. Finally, Fig.l shows that also for the ~ 30 targets 
in the visibility range of Arecibo (appearing at Tpp between 0.5 and 1.0 for the case of short or 
intermediate-duration pointings), SKAl-MID can go deeper when using 8-hr integrations (Arecibo 
can typically track sources for only 1—3 hrs). 

We note that a direct extrapolation from F pp to the expected number of pulsar discoveries in 
each GC is impossible for the majority of clusters, where there is currently no known pulsar. Such 
predictions are also uncertain for the GCs with known pulsars, because various effects (e.g. scin¬ 
tillation, fraction of binaries, Doppler smearing of the pulsar signal) can bias any simple relation 
between Tpp and the number of known pulsars. Nonetheless, despite all these caveats, the number 
of pulsars resulting from the direct multiplication of Tpp with the number of already known objects 
is compatible with the results of sophisticated statistical analyses (e.g. Bagchi et al. 2011; Chenna- 
mangalam et al. 2013) applied to some of the most populated GCs. These predict that the number 
of potentially observable pulsars should be (within 95% credible intervals) 82 — 259 in Terzan 5, 
48 — 137 in 47 Tuc, and 48 — 191 in M28. In summary, Tpp can be only taken as a figure-of-merit 
for the capabilities of SKAl-MID in the search for GC pulsars and is mostly useful for comparison 
with the current state-of-the-art Arecibo, GBT, and Parkes searches. In this respect, we note that a 
downgrade of SKAl-MID by a factor 50% in sensitivity would significantly reduce the effective¬ 
ness of the search for pulsars in GCs. In such a situation, it would be difficult to go significantly 
deeper than previous GBT searches, and impossible to go deeper than Arecibo. Thus, with a 50% 
reduction in SKAl-MID sensitivity, a significant sensitivity improvement would be limited to the 
GCs visible only from Parkes, i.e 1/5 of the total sample. 

3. Observing Pulsars in Globular Clusters with SKA2 

The great promise of SKA2 is illustrated in Fig. 2, where we show the distribution of the ‘com¬ 
pleteness percentages’ of a GC survey. For example, 50% completeness indicates the capability of 
probing half of the total area defined by the LF of a given GC — i.e. half of the area under the 
curve is probed. Even with 8-hr integrations, SKAl-MID cannot see all the potential pulsars in a 
given GC, whereas SKA2 may detect the majority or even all the active radio pulsars in a large 
sample of some tens of GCs — enabling detailed comparisons of their evolutions and structures, as 
well as an unbiased determination of the shape and cut-off of the LF of the millisecond pulsar pop¬ 
ulation. We also note that, when SKA2 will be operational, SKAl-MID will have likely detected 
at least one pulsar in most GCs, which will provide the possibility for SKA2 to speed up the search 

^Power-law LFs (with index between —1.0 and —0.7, and cut-off luminosity in the confidence range 0.1 — 
O.SmJykpc^) were also explored. Although the shape of the histogram distribution of Ypp is different with respect 
to the adoption of a log-normal function, the average value of Tpp over the whole GC population is similar. 
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SKA2-Mid - GC search completeness 



Figure 2: A histogram of the percentage completeness factors for GC pulsar searches with SKA2 (at 
1700 MHz), using 10-min, 2-hr, and 8-hr observations. See text for further details. 


for additional pulsars (and significantly reduce its computational cost) by using an already known 
approximate value of the DM. 


4. Conclusions 

Once constructed, SKA 1-MID and SKA 1-LOW will be the premier search machines for pul¬ 
sars in GCs. Their raw sensitivity will surpass that of the GET by a factor of about 3 and 4, 
respectively, which opens a large discovery space for faint sources outside of the Arecibo declina¬ 
tion range. Though the uncertainties are large, we estimate that a doubling or tripling of the known 
pulsar population is possible. In other words, anywhere from 100 — 300 pulsars can be found in 
a full census of ~ 150 GCs using 10-min integrations and only one tied-array beam synthesized 
from the SKA Core. This offers a great opportunity for early SKA pulsar science, even before all 
the collecting area is in place. A doubling or tripling of the known GC pulsar population means 
finding more exceptional pulsar systems for testing strong gravity, dense matter, and more. For 
example, only 3 of the known GC pulsars (all in Terzan 5) have spin rates greater than 500 Hz. A 
doubling or tripling of the known population gives good prospects for breaking the current 716-Hz 
rotation record held by Ter5ad. Besides the study of individual pulsar ‘jewels’ the ensemble of 
detected pulsars provide a unique probe of the dynamics and evolution. While SKAl can signif¬ 
icantly push such studies forward, SKA2 can bring us within a factor of a few of finding all the 
detectable (beamed towards us) pulsars in Galactic GC system, a population that has been estimated 
at 600 — 3700, and which is comparable to the total detectable Galactic field population. 
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